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Abstract 


A multispectral scanning spectrometer was used to obtain measurements of 
the reflection function and brightness temperature of clouds, sea ice, snow, and 
tundra surfaces at 50 discrete wavelengths between 0.47 and 14.0 Jim. These ob- 
servations were obtained from the NASA ER-2 aircraft as part of the FIRE Arctic 
Clouds Experiment, conducted over a 1600 x 500 km region of the north slope of 
Alaska and surrounding Beaufort and Chukchi Seas between 18 May and 6 June 
1998. Multispectral images of the reflection function and brightness temperature 
in 11 distinct bands of the MODIS Airborne Simulator (MAS) were used to de- 
rive a confidence in clear sky (or alternatively the probability of cloud), shadow, 
and heavy aerosol over five different ecosystems. Based on the results of indi- 
vidual tests run as part of the cloud mask, an algorithm was developed to esti- 
mate the phase of the clouds (water, ice, or undetermined phase). Finally, the 
cloud optical thickness and effective radius were derived for both water and ice 
clouds that were detected during one flight line on 4 June. 

This analysis shows that the cloud mask developed for operational use on 
MODIS, and tested using MAS data in Alaska, is quite capable of distinguishing 
clouds from bright sea ice surfaces during daytime conditions in the high Arctic. 
Results of individual tests, however, make it difficult to distinguish ice clouds 
over snow and sea ice surfaces, so additional tests were added to enhance the 
confidence in the thermodynamic phase of clouds over the Beaufort Sea. The 
cloud optical thickness and effective radius retrievals used 3 distinct bands of the 
MAS, with the newly developed 1.62 and 2.13 fim bands being used quite suc- 
cessfully over snow and sea ice surfaces. These results are contrasted with a 
MODIS-based algorithm that relies on spectral reflectance at 0.87 and 2.13 /im. 
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1. Introduction 

A knowledge of cloud radiative properties and their variation in space and 
time is especially crucial to the understanding of the radiative forcing of climate. 
High quality multispectral imagery acquired from high altitude aircraft or satel- 
lite platforms is the most efficient and reliable means of fulfilling these observa- 
tional requirements. Between 18 May and 6 June 1998, the NASA ER-2 high alti- 
tude research aircraft conducted eleven research flights over the north slope of 
Alaska and surrounding Beaufort and Chukchi Seas as part of the First ISCCP 
(International Satellite Cloud Climatology Project) Regional Experiment-Arctic 
Clouds Experiment (FIRE- ACE). 

The NASA ER-2 aircraft was equipped with seven sensors, the MODIS Air- 
borne Simulator (MAS; King et al. 1996), the Cloud Lidar System (CLS; Spinhirne 
et al. 1982), the Millimeter-wave Imaging Radiometer (MIR; Racette et al. 1996), 
the Airborne Multiangle Imaging Spectroradiometer (AirMISR; Diner et al. 1998), 
the Advanced Microwave Precipitation Radiometer (AMPR; Spencer et al. 1994), 
the High-resolution Interferometer Sounder (HIS; Smith et al. 1995), and the So- 
lar Spectral Flux Radiometer (SSFR; Pilewskie et al. 2000). The MAS was de- 
signed to obtain measurements that simulate those obtained from the Moderate 
Resolution Imaging Spectroradiometer (MODIS), a 36-band spectroradiometer 
launched aboard the Earth Observing System (EOS) Terra (King and Herring 
2000) and Aqua (Parkinson 2003) spacecraft. 

The strategy for FIRE-ACE included spaceborne remote sensing (polar or- 
biting satellites), high altitude remote sensing (NASA ER-2 at ~20 km), lower alti- 
tude remote sensing and in situ measurements (University of Washington CV- 
580, NCAR C-130Q, and Canada's AES CV-580 aircraft), ground-based meas- 
urements (radiation, clouds, meteorology, and surface fluxes), and modeling 
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studies (cf. Curry et al. 2000). FIRE-ACE took advantage of, and overlapped 
with, the Surface Heat Budget of the Arctic Ocean (SHEBA) experiment, which 
largely consisted of surface measurements on and around the Canadian Coast 
Guard icebreaker Des Groseilliers that drifted throughout the Beaufort and 
Chukchi Seas over a 12 month period. The ER-2 was based in Fairbanks and de- 
ployed over the north slope of Alaska and nearby Arctic Ocean and typically 
overflew the surface Atmospheric Radiation Measurement (ARM) site in Barrow, 
AK enroute to and from the SHEBA ice station. 

The main objectives of the ER-2 included: (i) comparing the spectral proper- 
ties of sea ice, tundra, and cloud layers, (ii) collecting MAS data to validate the 
MODIS cloud mask algorithm for distinguishing clouds from snow and sea ice 
surfaces in polar regions, (iii) collecting data for retrieving cloud radiative and 
microphysical properties over snow and sea ice surfaces during summer daytime 
conditions, (iv) determining radiative energy budget of clouds and sea ice in po- 
lar regions, and (v) comparing high-altitude remote sensing and ground-based 
observations of clouds and clear sky in polar regions, especially over the SHEBA 
and ARM long-term ground-based remote sensing sites. 

Clouds are generally characterized by higher reflectance and lower tempera- 
ture than the underlying surface. As such, simple visible and infrared window 
threshold approaches would appear to offer considerable skill in cloud detection. 
However, there are many surface conditions where this characterization of 
clouds is oversimplified and inappropriate, most notably over snow and sea ice 
surfaces. Additionally, some cloud types, such as thin cirrus and small cumulus 
are difficult to detect because of insufficient contrast with the surface radiance (or 
reflectance). Cloud edges cause further difficulty since the instrument field of 
view will not always be completely cloudy or clear. 

The 50 channel MAS offers an opportunity to explore and refine multispec- 
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tral approaches to cloud detection so that many of these concerns can be miti- 
gated, if not entirely eliminated. The intent of this paper is to describe the ap- 
proach and algorithms used to detect clouds during daytime conditions in the 
high Arctic during summer. This represents a subset of all conditions and chan- 
nels used to process global satellite data using MODIS. Given the results from 
the cloud mask, we have developed an algorithm, currently implemented in the 
MODIS global processing system, to estimate the thermodynamic phase of 
clouds (Platnick et al. 2003). Finally, we have retrieved cloud optical thickness 
and effective radius for the 'cloudy' scenes identified from the MAS during FIRE- 
ACE. Preliminary results obtained from the ER-2 on 4 June 1998 during FIRE- 
ACE are presented to illustrate the results of applying these cloud tests, thermo- 
dynamic phase decisions, and cloud microphysical retrievals to a wide range of 
conditions. 

2. Instrumentation 

The MAS is a cross-track scanning spectrometer that measures reflected solar 
and emitted thermal radiation in 50 narrowband channels. For the FIRE-ACE 
deployment, the configuration of the MAS contained channels between 0.47 and 
14.0 jum. Flown aboard the NASA ER-2 aircraft, the MAS scans through nadir in 
a plane perpendicular to the velocity vector of the aircraft (cross track), with the 
maximum scan angle extending 43° on either side of nadir (86° full swath aper- 
ture). At a nominal ER-2 altitude of 20 km, this yields a swath width of 37.2 km 
centered on the aircraft ground track, with a total of 716 Earth-viewing pixels ac- 
quired per scan. With each pixel having a 2.5 mrad instantaneous field of view, 
the spatial resolution is 50 m at nadir from the nominal aircraft altitude. 

Table 1 summarizes the band center and bandwidth characteristics as well as 
main purpose of each MAS band used to detect clouds and to retrieve their opti- 
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cal and microphysical properties during FIRE-ACE. Some of these bands are 
used to discriminate clouds from clear sky (cloud mask) where others are used to 
derive the optical and microphysical properties of clouds. The bands used for 
these purposes are identified in Table 1, and a detailed description of the algo- 
rithms will be presented in Section 3. Radiometric calibration of the shortwave 
(<2.5 jUm) channels was obtained by observing laboratory standard integrating 
sphere sources on the ground before and after this experiment, while calibration 
of the infrared channels was performed by viewing two onboard blackbody 
sources once every scan. A detailed description of the optical, mechanical, elec- 
tronics and data acquisition system design of the MAS can be found in King et al. 
(1996). 

3. Cloud Retrievals 

In order to derive cloud optical and microphysical properties, either from 
high-altitude remote sensing aircraft (NASA ER-2) or from satellite (Terra or 
Aqua), it is necessary to perform the following steps: (1) identify the probability 
of a given pixel being cloud contaminated, (2) determine the likely phase of the 
cloud to be analyzed, and (3) derive the cloud optical and microphysical proper- 
ties (such as optical thickness, effective radius, cloud top altitude, etc.). In this 
section we will describe each of these steps in turn, focusing on specific charac- 
teristics of the cloud retrieval process that are important considerations for 
clouds overlying snow and sea ice surfaces during summer, daytime conditions. 

a. Cloud mask 

The MAS cloud mask algorithm determines whether a given view of the 
earth surface is unobstructed by clouds or optically thick aerosol, and whether 
any clear scene is contaminated by shadow. The algorithm is divided into ten 
conceptual domains according to surface type (land, ocean, snow /ice, coastal, 
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and desert) and solar illumination (daytime and nighttime). In this paper, for 
application to MAS data over the high Arctic during northern summer, we are 
concerned only with daytime conditions, where daytime is defined for solar ze- 
nith angles 0o < 85°. The United States Geological Survey 1 km land / sea tag file 
is used for land /water discrimination. Over land, ecosystem classification is 
based on the 1 km International Geosphere-Biosphere Programme (IGBP) land 
cover dataset consisting of 17 land cover classes (Loveland and Bel ward 1997), 
and is supplemented by an automatic snow classification within the cloud mask 
algorithm for clear skies. Under conditions of cloudy skies and / or polar regions 
containing sea ice, this surface classification is supplemented by the National 
Snow and Ice Data Center (NSIDC) near-real time ice and snow extent (NISE) 
dataset that is available daily at 25 km resolution and is based on satellite-borne 
microwave radiometers. 

The threshold between a pixel being classified as cloudy or clear is some- 
times ambiguous due to instrument noise and natural variability. To allow for 
imprecise measurements of the real world and to accommodate a wide variety of 
applications, the cloud mask is more than a simple yes/ no decision. The cloud 
mask includes 4 levels of 'confidence' with regard to whether a pixel is classified 
as clear (cf. Ackerman et al. 1998; King et al. 1998, 2003). Confidence flags con- 
vey the strength of conviction in the outcome of the cloud mask algorithm tests 
for a given field-of-view (FOV). A confidence flag for each individual test, based 
on proximity to a threshold value, is assigned and used to derive a final quality 
flag determination for the pixel. The current scheme applies a linear interpola- 
tion between a low confidence clear (confidence = 0) and high confidence clear 
(confidence = 1) threshold for each spectral test. The final determination is a 
combination of the confidence indices of all applied tests. 

The confidence levels of each spectral test must be combined to determine 
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the final decision on whether a pixel is classified as clear or cloudy. Since several 
of the spectral tests are not independent of one another, the cloud mask algo- 
rithm first groups individual spectral tests. The MAS has many channels that are 
similar to or comparable to those in MODIS, but it lacks a few channels that are 
important in the MODIS cloud mask. On MAS, the band at 1.89 nm is in a strong 
water vapor absorption band, and is therefore quite comparable to the 1.38 nm 
band on MODIS. Both of these bands serve a similar purpose of aiding in the 
identification of subvisible, high-altitude cirrus clouds. MAS also lacks any 
bands in the 6.7 jtm water vapor region, a band that appears to be quite impor- 
tant in identifying water vs ice clouds (see below). Due to a lack of any band in 
the 6.7 jjm region, MAS performs its cloud mask spectral tests in four (rather 
than five) groups: (1) detecting thick high clouds with thermal infrared thresh- 
olds; (2) detecting thin high clouds with brightness temperature difference tests; 
(3) detecting low clouds using solar reflectance tests; and (4) detecting upper tro- 
pospheric thin clouds using near-infrared thresholds. Table 2 lists the MAS spec- 
tral tests that have been performed for all daytime conditions and for all five eco- 
systems. Thresholds for each group can be found in the MODIS cloud mask Al- 
gorithm Theoretical Basis Document (Ackerman et al. 2002). 

As described by Ackerman et al. (1998) and King et al. (1998), a minimum 
confidence is determined for each group as follows: 

Gj= i,n = m in[Fy]i=i,m, (1) 

where F\ is the confidence level of an individual spectral test, m is the number of 
tests in a given group, and N is the number of groups (4 in this case). The final 
cloud mask (Q) is then determined from the product of the results from each 


group. 
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Q = 





( 2 ) 


This approach is clear-sky conservative in the sense that if any test is highly con- 
fident that the scene is cloudy (F,- = 0), the final clear sky confidence is 0. The re- 
sult of this cloud mask test, together with processing path flags, heavy aerosol 
and shadow flags, and results from individual spectral and spatial tests, is stored 
in the 48 bit cloud mask file for each pixel (cf. Ackerman et al. 1998). In this pa- 
per these 2 bit cloud mask results will be interpreted as follows: (1) confident 
clear ( Q > 0.99); (2) probably clear (Q > 0.95); (3) probably cloudy (Q > 0.66); and 
(4) cloudy (Q <j 0.66). 


b. Thermodynamic Phase 

Once the cloud mask has been run for all pixels in a scene, it is necessary to 
interpret the results for application of the cloud retrieval code that follows (see 
below). Figure 1 shows the flowchart for determining the presence of clouds 
during the daytime using results reported in the cloud mask. Bit 0 of the cloud 
mask indicates whether a cloud mask result has been obtained for a given pixel, 
and if so, we check whether the pixel was observed during daytime or nighttime 
conditions. Next, we read the ecosystem classification bits in the cloud mask to 
determine which of the five ecosystems the pixel was acquired over. Since the 
ecosystem dataset applies to the static ecosystem of land vs water, supplemented 
by a dynamic snow index for clear sky pixels, it does not, at this point, incorpo- 
rate the dynamics of surface snow and sea ice that vary on a daily basis when the 
pixel is cloud covered. As this is of particular interest for our subsequent cloud 
retrievals, we supplemented this ecosystem map with the NISE ancillary dataset 
to determine whether snow or sea ice occurs beneath the potentially cloudy 
pixel. 
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Figure 2 shows the cloud mask tests that are run for pixels over snow and 
sea ice surfaces to determine whether the pixel contains clear sky, water cloud, 
ice cloud, or clouds of undetermined phase. If the scene is classified as cloudy or 
probably cloudy, then the tests available over snow and sea ice surfaces (cf. Table 
2) are run individually to estimate the thermodynamic phase of clouds in the 
scene. As there are limited tests available from the cloud mask, and the thresh- 
olds are set globally based on tests using MODIS data, it is possible that ambigu- 
ous phase identification still results from these cloud mask results. During the 
FIRE ACE experiment, all five ecosystems were encountered by the ER-2. Cloud 
mask tests for phase determination for each of these other ecosystems can be 
found at modis-atmos.gsfc.nasa.gov / MOD06JL2 / MAS_Cloud_Mask_Tests.pdf. 

Following the cloud mask tests, we run a number of additional tests aimed at 
increasing the confidence in the cloud thermodynamic phase and reducing the 
potentially large number of undetermined phase determinations (cf. Fig. 3). The 
first such test is a bispectral threshold algorithm based on the brightness tem- 
perature difference between 8.4 and 10.7 ,um. This is a simplified and robust 
modification of the trispectral phase algorithm first described by Strabala et al. 
(1994), which in turn is based on the difference in optical properties between wa- 
ter droplets and ice crystals, as described by Baum et al. (2000). This test is run 
only for those cloudy pixels for which the cloud mask tests led to an undeter- 
mined phase. Following this test all cloudy pixels are subjected to additional 
tests of the reflectance ratio between the shortwave bands at 2.1 and 1.6 jum and a 
visible, nonabsorbing band, at 0.66 Jim. The threshold values used over snow 
and sea ice surfaces are shown in Figure 3. For all other ecosystems, these 
threshold values are somewhat different and can be found at modis- 
atmos.gsfc.nasa.gov /MOD06_L2/ MAS_Final_Phase.pdf. 

In all significant ways, this is the same algorithm that is applied to determine 
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the thermodynamic phase of clouds using MODIS data from the Terra and Aqua 
spacecraft, except that the 1.6 pm shortwave infrared tests are not applied on the 
Aqua /MODIS instrument due to detector failures on this instrument. In addi- 
tion, for MODIS we implement additional cloud top temperature 'sanity checks' 
not implemented in MAS processing. These tests assure that any cloud whose 
cloud top temperature is warmer (colder) than 273 K (233 K) is liquid water (ice). 

c. Ice cloud and water cloud optical properties 

After the cloud mask and thermodynamic phase estimation has been per- 
formed, the physical and optical properties of ice clouds can, in principle, be re- 
trieved in a manner similar to that for water clouds, outlined below. Under the 
assumption of plane-parallel geometry, two other factors complicate the retriev- 
als of ice cloud properties (viz., the shape and orientation of the ice particles) that 
occur naturally in the atmosphere. Due to our limited knowledge accumulated 
thus far for ice clouds, the sensitivity of their retrieved properties on these two 
factors is still an ongoing area of current research. 

We have carried out single scattering calculations for ice crystals with differ- 
ent sizes and shapes. The sizes (measured along the maximum dimension of the 
crystals) range from 10 to 800 jum, while the shapes vary from bullet rosettes, 
hollow columns, plates, aggregates, to ice crystals with rough surfaces. The sin- 
gle scattering properties and phase functions at the 4 MAS bands for which cloud 
retrievals are performed (cf. Table 1) are presented for 12 size distributions ob- 
served for mid-latitude cirrus cloud systems. The single scattering calculations 
were carried out using the combination of the improved geometric optics 
method (GOM2) and the finite difference time domain (FDTD) technique (Yang 
and Liou 1996a,b). FDTD is applied to size parameter smaller than 20 whereas 
GOM2 is employed for size parameters larger than 20. 
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Following Takano and Liou (1989), we discretized the size distributions us- 
ing 5 size bins: 10-30 jum with bin center at 20 jum, 30-70 pm with bin center at 50 
pm, 70-170 Jim with bin center at 120 pm, 170-430 pm with bin center at 300 pm, 
and 430-1070 pm with bin center at 750 pm. The discretized number concentra- 
tion of ice crystals is depicted in Fig. 4 with respect to maximum dimension of 
the ice crystals. The effective size for a given size distribution is defined as fol- 
lows: 

D _ 3 fVn{D)dD _ 3 fWC/p,,, (3 ) 

6 2 J An(D)dD 2 AN ice ' 

where D, V, and A denote the maximum dimension, volume and projected-area 
of the ice crystals, respectively, A is the mean projected-area, N( ce is the total 
number concentration of ice crystals, IWC is the ice water content, and p , ce is the 
bulk density of ice. The preceding definition of effective size is applicable to ir- 
regular and complex ice crystal shapes, and does not require that the crystals be 
hexagonal columns, as in previous definitions. In addition, this definition of ef- 
fective size is directly related to IWC and crystal projected-area that can be meas- 
ured using a two-dimensional optical cloud probe. 

Wyser and Yang (1998) have shown that the effective size defined by Eq. (3) 
for a given size distribution is not sensitive to the detailed shape of the size dis- 
tribution. Finally, it can been shown that the present definition of effective size is 
a generalization of the effective radius defined by Hansen and Travis (1974) and 
given by 

00 M CO 

r e = J r 3 n(r)dr / J* r 2 n(r)dr, ^ 

0 0 

where n(r ) is the particle size distribution and r is the particle radius. This is the 
definition we have applied to our water cloud retrievals. In order to be consis- 
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tent, we have computed the effective radius of the 12 ice crystal models using Eq. 
(3) and dividing the resulting D e by 2. These values are contained within the in- 
dividual panels of Fig. 4. 

At this point, sufficient information concerning ice crystal habit and percent- 
age of individual habit for a given size distribution is not available. In this inves- 
tigation, as well as in processing of MODIS data, we assumed ice crystal habits 
and corresponding percentages for a model cirrus cloud based primarily on rep- 
licator data collected during FIRE-II-IFO. Figure 5 illustrates four shapes that we 
have numerically defined for single scattering calculations, where observed in 
situ ice crystal shapes are also presented for comparison. The percentage of each 
habit that we have assumed in our calculations varies with particle size. For 
small ice crystals (D <70 /urn), we assumed the particle shapes are relatively sim- 
ple, and characterized as 50% bullet rosettes, 25% hollow columns, and 25% 
plates. For large particles, we assume that bullet rosettes and aggregates domi- 
nate the particle size distribution, where the particles are composed of 30% ag- 
gregates, 30% bullet rosettes, 20% hollow columns, and 20% plates. The aspect 
ratio or dimensional relationship used in the calculations is taken from Auer and 
Veal (1970). For aggregates, the effect of crystal surface roughness is accounted 
for in the calculations, with the magnitude mathematically defined by using the 
Cox and Munk (1954) distribution. The variation from slight roughness to mod- 
erate roughness is randomly selected for each orientation by averaging the effect 
of particle orientations. For other crystal shapes, a smooth surface geometry is 
employed in the single-scattering calculations. 

In order to assess the impact of finite bandpass characteristics of the MAS on 
our cloud retrieval algorithm, we performed calculations of the phase function 
using 11-100 equally spaced wavelengths for each band. These results were then 
compared with computations based on a single wavelength at the bandpass cen- 
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ter. Numerical results show that the differences are negligible if the mean values 
of refractive index are used. For our computations, we used the complex refrac- 
tive indices of ice reported by Gosse et al. (1995) for wavelengths greater than 1.4 
jum, which deviate from data published by Warren (1984) by as much as 60% at 
some wavelengths. We use Warren's compilation for wavelengths below 1.4 pm. 

i 

For liquid water clouds, we have chosen to use the optical constants tabulated by 
Hale and Querry (1973) for bands below 0.872 fim, Palmer and Williams (1974) 

for the 1.618 |/m band, and Downing and Williams (1975) for the 2.133 fim and 

i 

greater bands. Table 3 shows the refractive index assumed in this paper for 
computations of both water and ice cloud properties. 

Figure 6 illustrates some examples of the phase functions of ice crystals at 3 
of the MAS wavelengths used in retrieving cloud optical properties during FIRE 
ACE. At 0.87 (im, a pronounced halo occurs at 22°, as expected. The magnitude 

| 

I of the forward peak decreases with increasing wavelength because the angular 

widths of the distribution of diffracted light is inversely proportional to size pa- 
rameter. The cirrus uncinus cloud (r e = 58.9 jum), shown on the right, has a more 
pronounced halo and forward peak than the cirrostratus model (r e = 14.45 jum) 

1 shown on the left. 

d. Retrieval of cloud optical thickness and effective radius 

The simultaneous retrieval of cloud optical thickness and effective radius is 
best achieved by simultaneously measuring the reflection function at a visible 
and near-infrared wavelength, and comparing the resulting measurements with 
theoretical calculations, as demonstrated by Nakajima and King (1990). This 
technique is especially accurate over dark ocean surfaces because the reflection 
function of the Earth-atmosphere system arises primarily from light scattering by 
the cloud layer, with little influence from the underlying surface. In comparing 
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measurements with theory, however, it is essential that the light scattering prop- 
erties of the cloud are modeled realistically, and that the cloud is properly as- 
cribed to either a water cloud or an ice cloud with corresponding optical proper- 
ties. For applications of this technique to polar regions involving clouds over 
snow and sea ice surfaces, it is necessary to further assume or otherwise estimate 
realistic values of the underlying surface reflectance at appropriate visible and 
near-infrared wavelengths. 

Figure 7 illustrates the underlying principles behind the simultaneous de- 
termination of x c and r e from reflected solar radiation measurements. The mini- 
mum value of the reflection function at each wavelength corresponds to the re- 
flection function of the underlying surface at that wavelength in the absence of 
an atmosphere. For the computations presented here, we assumed the sea ice 
present over the Chukchi and Beaufort Seas in May and June was Lambertian 
with a surface albedo of 0.60 at 0.87 ^m and 0.03 at 2.13 jra\, consistent with 
measurements obtained by the Cloud Absorption Radiometer during previous 
Arctic campaigns (Arnold et al. 2002). The dashed curves in Fig. 7 represent re- 
flection function contours for fixed cloud optical thickness, and the solid curves 
contours for fixed effective radius. These results show, as expected, that the re- 
flection function at 0.87 ju m is largely a function of cloud optical thickness, 
whereas the reflection function at 2.13 fim is largely sensitive to effective radius, 
with the largest reflectance values corresponding to the smallest particle sizes. 

As the arctic contains both water and ice clouds, we performed model cal- 
culations for both water clouds (Fig. 7a) and ice clouds (Fig. 7b) using the refrac- 
tive index, size distribution, and shape characteristics described in Table 3 and 
Section 3c. The data points superimposed on the theoretical curves of Fig. 7 cor- 
respond to MAS observations acquired at the observational solar and viewing 
directions specified in the figure for 4 June, a day in which the cloud phase algo- 
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rithm determined the scene in Fig. 7a as a water cloud and the scene in Fig. 7b as 
an ice cloud. 

Due to the relatively high surface albedo at 0.87 (im, we explored an alterna- 
tive technique for simultaneously retrieving cloud optical thickness and effective 
radius over snow and sea ice surfaces during FIRE ACE. This method, first de- 
scribed by Platnick et al. (2001), takes advantage of the fact that the surface al- 
bedo of snow and sea ice is quite low at 1.62 jum (cf. Arnold et al. 2002). Figure 8 
illustrates the theoretical relationship and corresponding data points of the si- 
multaneous retrieval of cloud optical thickness and effective radius using meas- 
urements at 2.13 jum and 1.62 jum, where we assumed the surface albedo of sea 
ice was 0.03 at both wavelengths. Although one loses the near-orthogonality of 
the retrieval of x c and r e using these two channels, and therefore one is more sen- 
sitive to calibration uncertainties in the measurements, the sensitivity to effective 
radius is much better because there is far less uncertainty in the value of the sur- 
face albedo. This technique is noticeably more robust for water clouds (Fig. 8a) 
than for ice clouds (Fig. 8b). 

4. Results from Observations 

During FIRE ACE, the ER-2 acquired 62.25 hours of data during eleven re- 
search flights between 18 May and 6 June 1998. These missions included re- 
peated measurements over ground-based radiometers, lidars, and radars at the 
SHEBA ice camp (8 flights) and ARM site in Barrow (11 flights) (Curry et al. 
2000). In what follows, we will describe results obtained from the ER-2 on one of 
these flights, demonstrating the performance of the cloud mask, thermodynamic 
phase algorithm, and subsequently cloud optical property retrievals over sea ice 
surfaces during the day. 
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a. Cloud mask and thermodynamic phase 

On 4 June, the ER-2 overflew the ARM site and then headed northwest to- 
wards the SHEBA ice camp to overfly cirrus being advected over the region from 
the northwest. It then flew four parallel north-south flight legs of approximately 
280 km in length, offset by 40 km, starting from an area east of the ice camp and 
ending over the station (76°51'N, 167°30'W). This flight was useful for cloud re- 
mote sensing and cloud mask validation, as well as sea ice retrievals in the pres- 
ence of changing cloud cover. Following the primary observations in the vicinity 
of the SHEBA ice camp, the ER-2 overflew the ARM site on its return leg over the 
Chukchi Sea, which was aligned with the ground track of the NOAA-14 polar 
orbiting satellite. The ER-2 pilot reported solid cloud cover throughout most of 
the mission, with clear sky observed only over the tundra between the Brooks 
Range and the Arctic Ocean coastline. Figure 9 illustrates the ER-2 ground track 
for this mission. 

The MAS acquired reflected solar and emitted thermal radiation measure- 
ments in 50 narrowband channels, 8 of which were used to test the cloud mask 
and determine the thermodynamic phase of clouds for various ecosystems, as 
outlined in Section 3 (cf. Table 1). Figure 10 shows images of Arctic stratus and 
altocumulus clouds 205 km in length over the Chukchi Sea some 120 km E of the 
SHEBA ice camp, where the ER-2 is flying from bottom (south) to top (north) up 
these images along the right-hand-most ground track parallel to the flight track 
over the SHEBA ice station. MAS visible and near-infrared channels can readily 
be used to detect clouds and sea ice, and to determine the thermodynamic phase 
and optical properties of these clouds. The panel on the left (0.87 pm) is sensitive 
both to scattering by clouds and to reflection by the underlying sea ice surface. 
In the lower portion of this image there are Arctic stratus clouds and a break, or 
hole, in the clouds. The center panel (1.62 pm) is largely sensitive to the thermo- 
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dynamic phase of the clouds (ice vs water) and is also useful for retrieving cloud 
optical thickness of water clouds over sea ice due to the low reflectance of sea ice 
at this wavelength (cf. Fig. 8). The water clouds that dominate the lower part of 
this image are somewhat brighter than the ice clouds that dominate the upper 
portion, and these ice clouds are themselves somewhat brighter than the under- 
lying sea ice surface. At 1.88 fim (right hand panel), the low-level water clouds 
are much darker than the mid-level ice clouds. Since this channel is in a strong 
water vapor absorption channel, these data are especially sensitive to cloud top 
altitude and are therefore of value in the cloud mask algorithm described in Sec- 
tion 3. 

These three channels give only a hint at the wealth of information available 
in 50 channel high spatial resolution data. Figure 11 shows a false-color com- 
posite image of this scene, together with images of the cloud mask, cloud mask 
tests, and final thermodynamic phase determination applied to this scene. The 
false-color image was constructed by contrast stretching and combining three 
spectral bands into one 24-bit image, where the spectral bands were assigned to 
red, green, and blue (RGB) 8-bit display channels. For this scene, the RGB as- 
signment was 1.62 jum (red), 0.87 jum (green), and 1.88 ^m (blue), corresponding 
to the same channels presented in Fig. 10. In this manner the mid-level ice 
clouds appear blue, low-level water clouds appear yellow or orange, and clear 
sky or sea ice appear green. 

The second panel in Fig. 11 represents the results of the cloud mask applied 
to this scene. As discussed in the previous section, the cloud mask assigns four 
levels of clear sky 'confidence' to each pixel. In this panel, the high confidence 
clear sky image on which the cloud mask is superimposed is the MAS reflectance 
in band 7 (0.87 jum). The cloud mask results are designated as cloudy, probably 
cloudy, probably clear, and high confidence clear. Comparing the cloud mask 
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results to the false color image on the left in Fig. 11, it appears that both the low- 
level water clouds in the lower portion of the image and the mid-level ice clouds 
in the upper portion of the image are being classified quite accurately. In the 
central portion of this image, on the other hand, the openings in the clouds are 
being correctly identified as clear sky while other portions of thin cloud around 
the openings are being classified as probably cloudy. 

Having determined the cloud mask for this scene, we ran the thermody- 
namic phase algorithm that makes use of individual cloud mask tests (cf. Fig. 2) 
to determine a first guess at the thermodynamic phase of the cloudy pixels. 
These results, shown in the third panel of Fig. 11, yield a high confidence of ice 
clouds in the northern half of the scene and liquid water clouds in the southern 
half of the scene, with a number of pixels of undetermined phase occurring pri- 
marily in regions where the cloud mask suggested probably cloudy pixels. 

Finally, we ran the bispectral threshold algorithm for undetermined phase 
pixels, and shortwave infrared tests for all pixels, as described in Section 3b and 
Fig. 3. The final result, shown in the right hand panel of Fig. 11, indicates that 
the upper portion of the cloud in this scene was identified as ice cloud while the 
lower portion of the cloud was identified as water, with very few remaining pix- 
els identified as undetermined phase. These phase results are used for the proc- 
essing path of all subsequent cloud optical property retrievals. 

b. Cloud optical and microphysical properties 

Having identified the corresponding scene as water or ice cloud, we per- 
formed cloud optical property retrievals as described in Section 3d. Figure 12 
shows retrievals of cloud optical thickness and effective radius derived using the 
retrieval algorithms illustrated in Figs. 7 and Fig. 8 for the solar and viewing ge- 
ometry appropriate to this scene, where the left hand pair of figures corresponds 
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to the 0.87 and 2.13 jum algorithm and the right hand pair of figures to the 1.62 
and 2.13 /urn algorithm. The data points superimposed on Figs. 7 and 8 corre- 
spond to southern and northern portions of the flight line that were identified as 
containing water or ice clouds in Fig. 11. As expected, the effective radius re- 
trievals yield rather similar results between the two algorithms, with the ice 
cloud particles being quite small in the arctic region, and where the smallest sizes 
generally occur for the optically thickest ice clouds in this scene. The failure of 
the 0.87/2.13 jitm retrieval in the optically thin northwest portion of this scene oc- 
curs because the measured reflection function at 0.87 /urn is less than the theoreti- 
cal calculations allow when Ag(0.87 gm) = 0.6. This is no doubt the consequences 
of leads and breaks in the sea ice in this region such that the actual underlying 
surface albedo is much less than 0.6. This only serves to emphasize the difficulty 
of using the traditional wavelengths to derive cloud optical thickness and effec- 
tive radius over a spatially nonuniform sea ice (or snow) surface. 

The cloud optical thickness in the traditional 0.87 and 2.13 ym algorithm is 
especially sensitive to assumptions on the underlying surface albedo, as illus- 
trated in Figure 7. King (1987) demonstrates that for optically thick clouds for 
which asymptotic theory applies (x c 2 : 9), the derived optical thickness depends 
explicitly on Ag such that uncertainties in surface albedo cause a systematic offset 
in derived optical thickness by an amount given by 4Ag/3(l-g)(l-Ag), where g is 
the asymmetry factor (-0.85 for water clouds). This offset is such that the re- 
trieved optical thickness retrieval for Ag = 0.5 (0.7) would be 4.44 greater (7.41 
smaller) than that for Ag = 0.6. Furthermore, as surface albedo increases to 0.7 
and beyond, an increasingly large number of pixels fail the cloud retrieval algo- 
rithm because the measured reflectance at 0.87 ym is less than the smallest values 
expected from theoretical calculations. As a consequence of uncertainties and 
variability in the underlying surface albedo of snow and sea ice surfaces at visi- 
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ble wavelengths, the Platnick et al. (2000) algorithm that combines two short- 
wave infrared bands (1.62 and 2.13 jim) is especially attractive for cloud optical 
property retrievals for water clouds. This technique is not, however, sufficiently 
sensitive to cloud optical properties of ice clouds due to the strong absorption 
characteristics of sea ice and ice particles. 

A numerical comparison of the two cloud retrieval algorithms can best be 
seen by examining histograms of retrieved cloud optical properties. Figure 13a 
shows comparisons of the marginal probability density function of cloud optical 
thickness for all water clouds in this scene, with Figure 13b showing the corre- 
sponding probability of effective radius. Figures 13c and d show corresponding 
probability density functions for all ice clouds in this scene. Though biases 
clearly exist in the cloud optical thickness retrieval, with the 1.62 and 2.13 fi m al- 
gorithm presumably the more accurate of the two, at least for water clouds, the 
effective radius is retrieved reasonably accurately in both cases. Based on this 
analysis as well as other retrievals using MAS data in the Arctic (not shown), the 
1.62 and 2.13 /im algorithm first described by Platnick et al. (2001) would appear 
to be a promising algorithm for routine cloud optical property retrievals over 
snow and sea ice surfaces during the daytime. This algorithm will soon be 
adapted to MODIS retrievals over these ecosystems for Terra data, though the 
1.64 jum band on Aqua /MODIS is not reliable. 

5. Summary and Conclusions 

High-resolution images of the spectral reflection function and thermal emis- 
sion of the earth-atmosphere system were obtained with the MODIS Airborne 
Simulator (MAS) operated from the NASA ER-2 aircraft during the intensive 
field component of the FIRE- ACE experiment, conducted over the Beaufort and 
Chukchi Seas of the Arctic Ocean between 18 May and 6 June 1998. Multispec- 
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tral images of the reflection function and brightness temperature at 10 wave- 
lengths between 0.66 and 13.98 fim were used to derive the probability of clear 
sky (or cloud), cloud thermodynamic phase, and the optical thickness and effec- 
tive radius of liquid water and ice clouds over sea ice in the high Arctic during 
summer. We compared two separate algorithms for determining the cloud opti- 
cal thickness and effective radius, one closely aligned to the version run opera- 
tionally to process MODIS data on the Terra and Aqua spacecraft, and the other 
a novel algorithm that shows promise, especially for liquid water clouds over 
snow and sea ice surfaces. 

We have shown that the new 1.62/2.13 fim algorithm is more robust at de- 
termining the cloud optical thickness and effective radius for water clouds over 
snow and sea ice surfaces, due primarily to the fact that both snow and sea ice 
have very low surface reflectance at these wavelengths. As a consequence, liquid 
water clouds provide a relatively strong reflectance contrast to the dark under- 
lying surface. This algorithm is less reliable for ice clouds due to the strong ab- 
sorption by ice clouds as well as the underlying snow and sea ice surface. 

In this paper we also describe the thermodynamic phase algorithm that is 
nearly identical to that implemented in global processing of MODIS data from 
the Terra and Aqua spacecraft. Both in MAS processing during FIRE ACE and in 
global processing of MODIS data, we characterize the optical properties of non- 
spherical ice clouds using 12 different size distributions, ice crystal habit that is a 
function of particle size, and state-of-the art light scattering models. These ice 
crystal characteristics, all described herein, have thus far proven to be robust in 
global processing of MODIS data and are likely to be extendable to more ad- 
vanced ice crystal models currently under development. The main difference at 
present is that MODIS utilizes the reflectance at 1.24 fim for retrieving cloud op- 
tical thickness over snow and sea ice surfaces, rather than 0.87 fim as described in 
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this paper. This wavelength, which is available on MODIS but not on MAS, has 
a smaller surface reflectance for dry snow (0.45) and wet snow or sea ice (0.25) 
than 0.87 jum, which is why it is chosen for MODIS processing of cloud optical 
and microphysical properties. In the future we plan to implement the 1.63 and 
2.13 Mm algorithm for Terra /MODIS processing of water clouds over snow and 
sea ice surfaces. 
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TABLE 1. Spectral and radiometric characteristics of all MAS bands used in the 
cloud mask and cloud optical property retrievals during FIRE ACE (day- 
time conditions). 


MAS 

Band 
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MODIS 

band 
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wavelength 

jum 

Spectral 
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lim 

Cloud 
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Cloud 

Retrievals 
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low thick cloud, shadow; cloud 
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low thick cloud, shadow; cloud 
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TABLE 2 . MAS cloud mask spectral tests executed (S) for a given processing 
path. 

Group 
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FIGURE LEGENDS 

Fig. 1. Flow chart for determining the presence of clouds during the daytime 
from the MAS cloud mask. 

Fig. 2. Cloud mask tests for an initial determination of the presence of water 
clouds, ice clouds, and clear sky over snow and sea ice surfaces during 
the daytime. 

Fig. 3. Flow chart for determining the final thermodynamic phase of clouds 
during the daytime over snow and sea ice surfaces. 

Fig. 4. Ice crystal size distributions for mid-latitude cirrus cloud systems. 

Fig. 5. (a) The ice crystal shapes commonly observed in cirrus clouds, (b) The 

microphysical model in terms of habit percentage that is used for 
MODIS/MAS retrievals of cirrus clouds. Note that D is the maximum 
dimension of an ice crystal. 

Fig. 6. Phase function computed for MAS bands 7 (0.87 /im), 10 (1.62 jum) and 
20 (2.13 jum) for the size distributions of Cs and Ci uncinus. 

Fig. 7 Theoretical relationship between the reflection function at 0.87 fim and 
2.13 fim for (a) water clouds and (b) ice clouds for various values of x c 
(dashed lines) and r e (solid lines) for specified values of surface albedo 
and 9o, 0, and (j). Data from measurements above arctic water and ice 
clouds are superimposed on the figure. 

Fig. 8 Theoretical relationship between the reflection function at 1.63 jum and 
2.13 fim for (a) water clouds and (b) ice clouds for various values of x c 
(dashed lines) and r e (solid lines) for specified values of surface albedo 
and 0o, 0, and <}>. Data from measurements above arctic water and ice 
clouds are superimposed on the figure. 

Fig. 9 Ground track of the ER-2 on 4 June 1998. 
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Fig. 10 MODIS airborne simulator (MAS) images of a 205 x 37 km section of a 
north-south flight line some 120 km due east of the SHEBA ice station 
on 4 June 1998 (5834 scan lines between 2047:00 and 2102:36 UTC). At 
0.87 jum it is difficult to distinguish the bright sea ice surface from 
overlying stratus or altocumulus clouds. At 1.62 fim, on the other hand, 
surface reflectance and thin ice clouds have lower reflectance than wa- 
ter clouds or thick ice clouds. At 1.88 fim, high-level cirrus clouds are 
readily distinguishable from the surface and low-level water clouds. 

Fig. 11. Composite image (left) of Arctic stratus, altocumulus, and sea ice on 4 
June 1998. The RGB assignment is red (1.62 Jim), green (0.87 /im) and 
blue (1.88 jum), and is based on the individual channels shown in Fig. 7. 
The second panel is the resultant cloud mask, the third panel the phase 
decision based on the cloud mask tests, and the final panel the final 
cloud phase decision after applying the bispectral and shortwave infra- 
red tests. The background image on which the cloud mask is overlaid 
is the MAS reflectance image at 0.87 ftm. The red lines in the false color 
image correspond to data points in the water and ice clouds overlaid on 
the theoretical curves in Figs. 7 and 8. 

Fig. 12. Cloud optical thickness and effective radius derived from the MAS im- 
ages on 4 June 1998. The pair of images on the left correspond to the 
0.87 and 2.13 jum algorithm, and the pair of images on the right corre- 
spond to the same scene analyzed using the 1.62 and 2.13 jum algorithm. 
In all cases, we have used a split color bar to designate those pixels 
analyzed using the liquid water and ice cloud processing path. 

Fig. 13. Marginal probability density function of cloud optical thickness and ef- 
fective radius for all water (a and b) and ice (c and d) pixels in the MAS 
flight line on 4 June 1998. The pair of probability density functions in 
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each panel correspond to the probability distribution of cloud retrievals 
for the 0.87 and 2.13 jum algorithm (dotted distribution) and the 1.62 
and 2.13 jum algorithm (shaded distribution). The distribution arising 
from changing the surface albedo at 0.87 jim from 0.6 to 0.65 is shown 
as a solid curve in all panels. 
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See Fig. 2 


Figure 1. Flow chart for determining the presence of clouds during the daytime 
from the MAS cloud mask. 
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Figure 2. Cloud mask tests for an initial determination of the presence of water 
clouds, ice clouds, and clear sky over snow and sea ice surfaces dur- 
ing the daytime. 
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Figure 4. Ice crystal size distributions for mid-latitude cirrus cloud systems. 
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Figure 5. (a) The ice crystal shapes commonly observed in cirrus clouds, (b) 
The microphysical model in terms of habit percentage that is used for 
MODIS/MAS retrievals of cirrus clouds. Note that D is the maximum 
dimension of an ice crystal. 
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Figure 6. Phase function computed for MAS bands 7 (0.87 jum), 10 (1.62 ium) and 
20 (2.13 ium) for the size distributions of Cs and Ci uncinus. 
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Reflectance (0.86 pm) 

Figure 7. Theoretical relationship between the reflection function at 0.87 jum 
and 2.13 for (a) water clouds and (b) ice clouds for various values 
of x c (dashed lines) and r e (solid lines) for specified values of surface 
albedo and 0o, 0, and 4>. Data from measurements above arctic water 
and ice clouds are superimposed on the figure. 
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Reflectance ( 1 .62 jim) 

Figure 8. Theoretical relationship between the reflection function at 1.63 jum 
and 2.13 i«m for (a) water clouds and (b) ice clouds for various values 
of x c (dashed lines) and r e (solid lines) for specified values of surface 
albedo and 0o, 8, and <j). Data from measurements above arctic water 
and ice clouds are superimposed on the figure. 
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Figure 10. MODIS airborne simulator (MAS) images of a 205 x 37 km section of a 
north-south flight line some 120 km due east of the SHEBA ice station 
on 4 June 1998 (5834 scan lines between 2047:00 and 2102:36 UTC). At 
0.87 jum it is difficult to distinguish the bright sea ice surface from 
overlying stratus or altocumulus clouds. At 1.62 jum, on the other 
hand, surface reflectance and thin ice clouds have lower reflectance 
than water clouds or thick ice clouds. At 1.88 fim, high-level cirrus 
clouds are readily distinguishable from the surface and low-level wa- 
ter clouds. 
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Figure 11. Composite image (left) of Arctic stratus, altocumulus, and sea ice on 4 
June 1998. The RGB assignment is red (1.62 jum), green (0.87 jum) and 
blue (1.88 jum), and is based on the individual channels shown in Fig. 
7. The second panel is the resultant cloud mask, the third panel the 
phase decision based on the cloud mask tests, and the final panel the 
final cloud phase decision after applying the bispectral and shortwave 
infrared tests. The background image on which the cloud mask is 
overlaid is the MAS reflectance image at 0.87 fim. The red lines in the 
false color image correspond to data points in the water and ice 
clouds overlaid on the theoretical curves in Figs. 7 and 8. 
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Figure 12. Cloud optical thickness and effective radius derived from the MAS 
images on 4 June 1998. The pair of images on the left correspond to 
the 0.87 and 2.13 jum algorithm, and the pair of images on the right 
correspond to the same scene analyzed using the 1.62 and 2.13 |Um al- 
gorithm. In all cases, we have used a split color bar to designate those 
pixels analyzed using the liquid water and ice cloud processing path. 



Probability Density Function Probability Density Function 


KING ET AL. 


44 



Cloud Optical Thickness Cloud Effective Radius (pm) 


Figure 13. Marginal probability density function of cloud optical thickness and 
effective radius for all water (a and b) and ice (c and d) pixels in the 
MAS flight line on 4 June 1998. The pair of probability density func- 
tions in each panel correspond to the probability distribution of cloud 
retrievals for the 0.87 and 2.13 fim algorithm (dotted distribution) and 
the 1.62 and 2.13 algorithm (shaded distribution). The distribu- 
tion arising from changing the surface albedo at 0.87 fim from 0.6 to 
0.65 is shown as a solid curve in all panels. 






